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SUMMARY 


On the basis of linear theory, the zero-lift wave drag of a particu- 
lar family of imswept, tapered wings with linearly varying thickness ratio 
and symmetrical parabolic-arc sections has been calctilated. The case of 
the wing with a given root, thickness ratio is given primary consideration 
in this paper with the view toward its use for missiles with all-movable 
fins where the root thickness must be large enough to allow for a rigid 
attachment to the trunnion and controlling mechanism. By con^jaring the 
drag for these wings with that for a corresponding constant-thickness- 
ratio wing with rhombic sections, it is found that the variable-thickness- 
ratio wings can be used to advantage with no serious structiiral penalties 
if the wings are assimied to have the same given TOOt thickness ratio or 
the same internal volume. 


INTRODUCTIOJ 


Zero-lift drag calciilatlons have generally not been made for tapered 
wings with c\rrved surfaces because the thickness functions, from which 
the source distributions are obtained, are us ua lly of such a nature that 
the drag equations are nonlntegrable, although in some cases pressure 
distributions may be foimd. References 1 and 2 present a numerical method 
for computing the pressure wave drag of delta and arrow plan-form wingsj 
reference 1 is for biconvex, constant-thlckness-ratlo sections, and. ref- 
erence 2 determines the TniniTmirn wave drag for constant-thickness-ratio 
wings. 

By modifying the equation which describes the surface of an unswept, 
tapered wing of biconvex section and constant thickness ratio, for which 
the drag has not yet been found analytically, an unswept, tapered wing of 
biconvex section and linearly varying thickness ratio is obtained for 
which the drag can be found. This is not a coopletely general wing in 
that the taper ratio is equal to the ratio of the tip thickness ratio to 
root thickness ratio. Nevertheless, it is of immediate practical Interest, 
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particularly as applied to missiles with all-movahle controls for which 
the primary consideration often is that the root thickness be large 
enough to allow for a rigid attachment of the wing to the trunnion and 
controlling mechanism witho-ut incurring any large stress concentrations. 
Since missile wings are often subjected to large, instantaneous deflec- 
tions, the Jxmcture of the wing ajid trunnion must be well designed. 

The iinswept wing was chosen for this investigation because, for the 
hi^er Mach numbers at which missiles and future aircraft will be oper- 
ated, the unswept wing is generally sr^jerior to the swept wing as far as 
drag is concerned. As an Illustration, figure 1, vhlch is taken from 
figure 10(b) of reference 3 , is presented. Figure 1 presents the ratio 
of the drag coefficient of a rectangular wing to the drag coefficient 
of an xmtapered, sweptback wing of the same aspect ratio (A = 5), plan- 
form area, thickness ratio, and parabolic-arc section as a function of 
Mach ntmiber for varloois angles of sweepback. After the leading edge of 
the sweptback wing becomes supersonic, the drag characteristics of the 
rectangular wing are seen to be as good as or superior to those of the 
sweptback wing. As the Mach nimtber increases, the angle of sweepback 
necessary to reduce the drag below that of the rectangular wing becomes 
structurally unfeasible; therefore, the rectangular wing is preferable 
to the sweptback wing at the hi^er MaCh numbers. 

For a given thickness at the root, which is the primary conceni in 
this paper, an unswept, tapered wing will, have better drag character- 
istics than a rectangular wing of the same aspect ratio and plan-form 
area because the tapered wing will have a larger root chord and, conse- 
quently, a lower thickness ratio. A further decrease in the drag can be 
obtained by retaining the tapered plan form and the given root thickness 
ratio but reducing the local spanwise thickness ratio. Of course, this 
spanwlse reduction in thickness ratio Increases the local bending stresses 
so that, for any particular application, a satisfactory compromise between 
allowable stresses and drag reduction miist l)e reached. 

Since the higher Mach ninnbers are of interest and the wing of the 
present paper is \mswept, the drag calculations have been made only for 
the case of supersonic leading edges. 


SYMBOLS 


A 


aspect ratio. 


ke^/s 


B = PA 


C 


D 


wave-drag coefficient, D/qS 
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Cjjjj wave-drag coefficient of rectangular wing 

Cjjg wave-drag coefficient of sweptback wing 

c(y) local chord 

Cj, root chord 


D 

F 

1 - 
M 

m(y) 

N 


Ap 


wave drag 
frontal area 

moment of Inertia about local chord 


20 

Mach number 

local bending moment 

c^(l - A) 

sweepback parameter, , tangent of leading-edge 

sweepback angle 

local pressure minus free-stream pressure 

1 p 

dynamic pressure, 



S 

S' 


s 

t(y) 


plan-form area 
area of integration 
semi span 

maximum local thickness of wing 
maximum root thickness 


skin thickness 
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V free-stream. velocity 

I 

Vj_ internal volume 

V vertical perturbation velocity 

x^y^z Cartesian coordinates 

p = - 1 


A 


P 

c^(y) 

T(y) 


leading-edge sweepback angle 


. . , . Tip chord 

wing taper ratio, — 

^ ^ Root chord 

coordinates of source point 

free-stream density 

maximimi local bending stress. 


t(y) 

local thickness ratio, - 7 — r 

c(y) 


^xx(y) 


Tj, root thickness ratio, — 

'^r 

0 perturbation velocity potential 

Unless otherwise denoted, primed symbols refer to the constant- 
thickness-ratio wing. 


ANALYSIS 


The siirface of the variable -thickness -ratio wing of the 
investigation has parabolic-arc sections as described by the 


present 

equation 


z 


= 2Tj.Cj. 




( 1 ) 
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The maximum-thickness line is at x = Cp/2. The projection of the 
maximum-thickness line on the yz-plane is descrihed hy 


z 


Ir^r A 



and the thickness ratio is 



By using the definition of N, equation (2) may be rewritten as 

T(y) = T^[l - (1 - 

from which 



= A 


( 2 ) 


(5) 


Thus, the wing which is treated herein is not completely general in 
that the ratio of the tip and root thickness ratios is directly related 
to the taper ratio. Nevertheless, equation (l) does allow the drag to 
be calculated for cvirved-s\arface airfoils of practical interest, as men- 
tioned previously. 


The linearized partial-differential equation for the perturbation 
potential in steady supersonic flow is 


2 8^0 8^0 ^^0 
Sy^ ^z^ 


(5) 


If the disturbances are assumed to be small, the boundary conditions 
on the surface of a thin, synmetrical wing at zero angle of attack may be 
satisfied to the first order in the plane of the wing. Hence 



( 6 ) 


Reference 4 shows that -a solution of eqoiation (5) which satisfies condi- 
tion (6) is, for a source distribution. 
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v(l,T}) d| dT) 


(7) 


where w(|,tj) = 


az(l,Ti) 

d| 


unit area. 

To the- first order. 


■V is proportional to the source strength per 


the pressure in the plane of the wing is given 


hy 


The drag is expressed as 


<1 vW4=o 



Ap dz 
q dx 


dx dy 


( 8 ) 


(9) 


and can be found by either of two essentially equivalent approaches. 

The first method follows from reference 5 and is the more general 
of the two procedures. Since, for the zero-lift case, the potential in 
the z = 0 plane off the wing is zero, the potential at any point on 
the wing is the integrated effect on that point of all the sources within 
the wing boundaries which can influence that point. Thus, the wing of 
figure 2 has five regions of influence. The potential for each region 
can be found by integrating equation (j) over the area bounded by the 
wing plan form and the Mach lines which are drawn forward from the 
arbitrary point of each region. In practice, however, it is necessary 
to perform the integration only for the arbitrary point of region (5)^ 
since the potential for each of the other four regions consists of the 
appropriate real parts of the potential of region ( 5 ) . The drag can 
therefore be found from equation (9) t>y integrating the product of the 
slope and the pressure for region (5) over the whole plan form and taking 
the appropriate real parts of the result. 

The second method consists of the superposition of source distribu- 
tions, a procedure which is described in reference 6. Although this 
method is not so straightforward as the first, its use is often advanta- 
geous, particularly if some previously calculated results are applicable, 
as was the case with the present paper. As shown in figure 5; there are 
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three regions over which equation ( 9 ) is 
cases to the condition that B > K, that 
leading edges). These regions are: 


applicable, subject in adl three 

is pA > ^ (supersonic 

1 + X 




Region I; B ^ 2 

Region II: — ^ ^ B < 2 

1 + A 

Region III: 0 < B < — ^ 

" " 1 + A 

d 


B > K 


Thus, if a wing has a configuration with respect to the Mach lines that 
fits region II, the fxonctions of equation ( 9 ) relating to regions I and II 
are real, and those for region III are imaginary. 

The drag functions obtained from equation ( 9 ) have been evaluated by 
the second method and are given in the appendix for the three different 
regions. The drag is plotted in figure 4 as a fxmction of B = pA and A. 


RESUITS AUD DISCUSSIOH 


In addition to having good drag characteristics, a practical wing 
must be structurally sound. In order to ascertain the relative drag and 
struct-ural properties of the present wing, the drag and maximum local 
bending stresses of the variable-thickness-ratio wing are compared with 
the same properties of a corresponding constant-thlckness-ratio wing. 

The results for the constant -thiclsness-ratio wing which has rhombic sec- 
tions are obtained from reference J. 


Given Root Thickness Ratio 

Since the main application of the present wing is to missiles with 
all-movable fins where the thickness at the root must be sufficient for 
a rigid attachment of the wing to the trunnion and controlling mechanism, 
the drag and stiructural characteristics of a variable- and a constant - 
thickness-ratio wing of the same plan form should be made on the basis 
that they both have the same root thickness ratio. If primed symbols 
denote values pertaining to the constant-thickness -ratio wing (obtained 
frcm ref. T) ajcid unprimed symbols represent the variable-thlckness-ratio 
wing, the ratio of the drags of the two wings on the basis of a given 
root thickness ratio is 
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Cp CpPAr^ 

Cp’ CD«3/Tr‘^ 


( 10 ) 


Equation (lO) is presented in flgiare 5^ from which is seen to 

be a strong function of 'X, with drag reductions obtained for values of 
X less than about 0.70. 


It shotild be noted that the constant-thlckness-ratlo wing used for 
comparison has rhombic sections, which makes Cjj^Crp’ less favorable 

than if it had parabolic-arc sections (that is, for a given thickness 
ratio, a rhombic section has less wave drag th^ a parabolic-arc section) . 
According to reference 8, if a constant-thickness-ratio wing has 
parabolic-arc sections, all values of Cjj/Cjj' should be snultiplled by 


0 . 75 . 


As indicated by linear theory, two identical plan forms at the same 
angle of attack will have the same lift distribution and, consequently, 
the same local bending moment. Since the may j mum local bending stresses 
are 


. 1 2khM (11) 

the ratio of the maximum local bending stresses of the solid variable- 
thlckness-ratio wing to the constant-thlckness-ratio wing is 


Now, 


and 


(y(y) ^ •fa(y)l XX* (y) 

c^*(y) ixx(y)'t'(y) 


t(y) = T(y)c(y) 


I 


XX 


oT.E.(y) Pz(x,y) 

2 / / z^z dx 

OL.E.(y) Oq 


( 12 ) 


( 15 ) 


( 1 ^) 


Therefore, 
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g(y) 

tJ'(y) 



which, on the basis of given root thickness ratios, becomes 


(15) 


g(y) 

5'^y^Jsolid 



(16) 


Implicit in equation (15) is the assumption of solid wings. For 
large wings where, in order to reduce the weight, a hollow wing construc- 
tion is used, calculation of based on the assumption that the skin 

carries all the load, will probably more nearly represent the true state 
of affairs than does equation (l5). For a wing with constant skin thick- 
ness tsj where is the point on the x-axls where hollowness begins. 


^x-,(y) pz(x,y) 


c/2 pz(x,y) 


IjQC / z az dx + 4 / / z^dz dx 

'-/L.E.(y) -Jo ^X 3 _(y) '-^z-(x,y)-tp 


(17) 


frcm which, by using equation (12), 


g(y) 
cr‘ (y). 


55rT^ 


hollow 8 t/[i - (1 - A)|] 


1 - "" 


1 . r 

% 

1 - (1 - 

. 

1 - (1 - A)| 

[i - (1 - x)|]^ 


1 - 


1 - 


2r 


1 - (1 - 


W 


( 18 ) 

where r = — ^ 

Tr - 2 
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If the term to the 7/2 power is expanded and r Is a3JLowed to 
approach zero, then 


g(y) 

, 2f , ,v' 

56 t ^ 1 - (1 - >.)| 


(19) 


which, for given root thickness ratios, is 


g(y) 

ff'(y) 

hollow 


39 

56 [1 - (1 - >v)J 


( 20 ) 


Eqmtions (16) and (20) are plotted in. figure 6, where the solid lines 
are for the solid wing and the dashed lines are for the hollow wing which 
has a skin thickness that approaches zero. The hatched portion between 
each set of solid and dashed lines represents varying .degrees of skin 
thickness. From this figure it may he seen that a hollow variahle- 
thickness-ratio wing generally ccMpares more favorably with a hollow 
constant-thickness-ratio wing than the solid wings do. 


A comparison of figures 5 and 8 reveals that, as the taper ratio is 
decreased below about 0.70, increasingly large drag reductions are 
obtained but with correspondingly large increases in the ratios of local 
bending stresses, especially in the vicinity of the tip. Within limits, 
however, this Increase in outboard bending stresses is not serious. In 
figure 7 ratio of local bending stresses to root bending stresses of 
a solid wing with a constant thickness ratio is plotted against the span- 
wise position, as an illustration. The curves are based on the conserva- 
tive assumption of constant pressure distribution. It may be seen that, 
for all but the case of A = 0, the outboard sections are hl^ly under- 
stressed. That is, from a structural standpoint, in all but the case of 
A = 0, a constant -thickness-ratio wing wastes material. 


As a matter of interest, the drag and bending moment of the two 
wings can also be compared on the basis of given frontaJL area and given 
volume. 


Given Frontal Area 

The projected frontal area of a wing is 



( 21 ) 
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from, which 


F = 




1 + A + 



( 22 ) 


F' 


Tj.'Cj,8 

“TT“ 


(1 + A) 


(23) 


For F = F' with all plan-form dimensions the same, 

^ 3xr« (1 + A) 

^ 2(1 + A + A^) 

Substituting equation (24) into equations (lO), (l5)^ (3-9) ^ on 

the basis of a given frontal area, yields 



% « 9(1 + A)^ 

*^d' 4(1 + a + A2)2 Qj^TpyVj.»2 


g(y) 
o’ (y) 


-* solid 


35(1 + A + a2)‘ 


144(1 + A)^[i - (1 - 


(25) 


( 26 ) 


o(y) 
o' (y). 

hollow 


35 (1 + A + 



126(1 + A)^|i - (1 - A)fj 


(27) 


Equations (25), ( 26 ), and (27) are plotted in figures 8 and 9 axid show 
that, for all taper ratios, the drag of the variable -thickness -ratio wing 
is higher than the drag of the constant-thickness-ratio wing, whereas the 
stresses in the wing are lowered. 
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Given Volume 

The internal volume of a wing is 

>s /-.T.E.(x) 


from which, 


n T.JEi. 

z(x,y)dx dy 

..E.(x) 


Vi = (i + 7^ + 7,2 + ^5) 


(28) 


(29) 


Vi« = + X + A^) 


(50) 


For Vi = Vi' with all plan-form dimensions the same. 


Tr = 


Ty'(l + A + A^) 
1 + A + A^ + 


(51) 


Substituting equation (51) into equations (lO), (l5)^ and (19); on 
the basis of a given internal volume, yields 


Cp _ (i + A + A^f CppAr^ 

’ (l + A + a 2 + A^f Cp‘p/x/2 


(52) 


ra(y) 1 

(y)J 


55(1 + A + A^ + A^)' 


solid 64(1 + X + a2)^ ji - (1 - X)f ^ 


(55) 


ra(y) 1 

LcT’(y)J 


hollow 


55(1 + A + a2 + X^)^ 

56(1 + A + A2)^ji - (1 - A)| 


(5i^) 
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Equations (52), (55) ^ and. (3^) are plotted in figures 10 and 11. 

These figures show that drag reductions are obtained for all taper ratios, 
whereas the structural characteristics are somewhat worsened by comparison 
with the constant-thickness-ratio wing. Reference to figure 7^ however, 
shows that, within limits , these poorer structural characteristics are 
tolerable. 


CONCinDING REMARKS 


On the basis of linear theory, the zero-lift wave drag of an unswept, 
tapered wing with linearly varying thickness ratio and parabolic-arc 
sections has been calculated. Although this wing is not completely 
general, in that the ratio of the tip thickness ratio to root thickness 
ratio is equal to the taper ratio, it is nevertheless of immediate practi- 
cal interest for such applications as missiles with allnmovable fins where 
the primary consideration is often that the thickness at the root be sirC- 
ficient to allow a rigid attachment of the wing to the trunnion and con- 
trol mechanism. 

The drag and bending-stress characteristics have been ccmpared with 
the corresponding characteristics of a constant-thickness-ratio wing of 
the same plan form but with rhombic sections bn the basis of the same 
root thickness ratio, the same frontal area, and the same internal volxmie. 

For the case of the same root thickness ratio, the primary concern 
of this paper, the ratio of the drag of the variable-thickness-ratio wing 
to the drag of the constant-thickness-ratio wing becomes less than 1.00 
at a taper ratio of about O. 7 O. Decreasing the taper ratio decreases 
the drag ratio but increases the ratio of local bending stresses of the 
variable-thickness-ratio wing to the local bending stresses of the 
constant-thickness-ratio wing. This bending-stress ratio is a maximuTn 
at the tips. However, inasmuch as most constant-thickness-ratio wings, 
when designed for a given root bending moment, will be under stressed 
oxjtboard, an outboard increase in bending stress is not too serious, 
within limits. 

Although the wings with the same frontal area were better structurally 
than those in the preceding case, the drag was adversely affected for all 
taper ratios. 

For the wings with the same internal, volimoe, the drag ratio is less 
than 1.00 for all taper ratios and, although the outboard stresses are 
increased sli^tly over those for the case of the same root thickness 
ratio, a limited Increase in outboard bending stresses is not 
objectionable . 
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It appears, therefore, that for the cases of the same root thickness 
ratio and the same internal volume, a satisfactory ccaaprcmiise between drag 
and stress considerations can be reached which will allow the present wing 
to be used to advantage. 


Langley Aeronautical Laboratory, 

National Advisory Cormnlttee for Aeronautics, 
langley Field, Va., January IJ, 1955* 
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APPENDIX 


DRAG FUNCTIONS 


The drag fimctions (eq. ( 9 )) have been evalmted for the three regions 
indicated in figure 5 where pA > ^ or B > K. These functions 


1 + X 


are; 


For B ^ 2, 


/CdP\ 128 b 

Fb^ - 2B^^ + lOB^K^ - 1)-K^ 

VTr^/j A)^\/b2 - k2 

[ 12KI 

[b2 _ Y?) 

? 


- 2 K^) 
(b 2 - K^) 


12K' 


-1 K ^ - 28B^K^ + lOK^ 

cos ~ + =-T^ 

® 56 (b 2 - 


ixaHe^ - 2K^) A^ 
12 k(b^ - K^) 12^fe2 - 


For 


2A 


1 + A 


^ E < 2, 


(Al) 


VtViI Vr /i ^ ;^j2 _ 


K® 


I 28 B - ac^) B^(i + \) - i»(i - _ 


12 ?(b^ - K®) 




(1 - 7.)5\/b2 -K^ ^ 

2 + \Ik -3^ 

+ ? 

6k2 ^ 

B 

. 


B^(tB^ - 22B^ + 


12 K(b® - JS^'f 


s\(.SB^ - K^) ^ - 1^) _ 2^ 

3 k(b^ - K^) 2 k(b^ - K®) 


-1 B^(1 + X) + 4(1 - 7^) , 


4£ 


-14B^ + 6tb^I^ - 25k‘*' . ^ ?>(2SB^ - 23K^) _ ITX^ 
56 (b^ - 56 (b^ - K?)^^ 



_ 

•* 

ITX^ 

7.5 

(1 + \)\lk - b^ 

36 

12\/b^ - 

4 


(A2) 
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For 


0 ^ B < ■ 

1 + A 


> 


f 


Jr /III Vr/ii *(1 + 


- 8 B^ -I- ^ IgB^K**' + ^ 

12 K(b 2 - k 2 )^ 


2X(a2 - K^) _ Tv^CsB^ - K^) ^ ^ b2(1 ^ + 4X(1 - X) ^ 

5k(b2 - K^) 2K(b2 - B?) 3K 3(1 + X) 


ICB^ - 37BgKg + izsd' _ xC^B*'’ + kOB^ - ^ 

36(b2 - r2)^ ^ 72K?(b^ - K?)^ 


X^(5B^ + 5K^) 
18 K? - K? 


7? 

ikT? - B^(1 + X)® 

12K® 



^(a + x^) - b 2 (i - x)^J iQgg 


2X 


d 


Ki? - B^(l + X) 
B(1 + X) 


These equations are indeterminate for the case of A = 1. In 
limit as A — >1, they reduce to 

VrVr 5 




52 

35T 


sin 


- B A " ^ b(6 - B^) 


log 


1 + \/iT 


B" 


B 


(A3) 

the 

(A5) 

(a6) 
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Flgiire 1 .- Variation of ratios of wave-drag coefficient of rectangular 
wings to wave-drag coefficient of imtapered, sweptback wings plotted 
against Mach number. Both wings have same plan fonn, parabolic-arc 
sections, and aspect ratio (a = 5 )* 



i 



Figure 2.- Regions of Influence for determining presBure distribution of 

an tmswept, tapered wing. 
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figure 4.- Veurlatlon of zero-lift wave-drag coefficient of tmavrept wings 
having linearly varying thickness ratios and parahollc-arc sections 
with pA for various values of taper ratio. 
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•atlo. 
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Figure 8.- Variation of ratios of vave-drag coejfflcients of varla'blfe- 
tMctaie SB -ratio wing and constant-thickreBB-ratio wing plotted against 
j3A for various values of taper ratio. Both wings have same frontal 
eirea. 












figure 10.- Variation of ratio of vrave-drag coefficients of variable- 
thlcknesB-ratio wing and conetant-thickness-ratio wing plotted against 
pA for various values of taper latio. Both wings have eajne internal 
volume. 






